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Abstract : This paper examines the ability of a 

real-time “disturbance-observer” to adapt-to and 
closely estimate the time-behavior of a 
disturbance-input w(t) (and of it’s state-vector z(t)) 
even when t he actual w(t ) time-behavior deviates 
from the observer’s “internal-copy” of the 
nominal/predicted H>(f)-behavior. By means of 
technical explanations and confirming simulation 
studies of numerical examples, the disturbance- 
observer’s adaptive ability is explained in terms-of 
the underlying spline-model used to derive the 
disturbance state-model and the intrinsic dynamic 
characteristics of a state-observer. 

I. Introduction 

The modern-control concept of a “disturbance- 
observer” for linear dynamical systems was 
originally developed in a 1966-68 research-study 
funded by NASA’s Marshall Space-Flight Center 
through the Huntsville Office of the General 
Dynamics Corporation [1], [2]. Those earliest 
results, obtained by variations of the optimal Linear 
Quadratic Regulator control-design methodology and 
by a novel alternative, purely (linear) algebraic/ 
stabilization control-design method [2; pp. ], were 
first published in a series of 1968-1970 journal 
papers [3]-[5], 

The concepts of: (i) a disturbance-“state” z(t ), (ii) 
an associated disturbance state-model z = Dz 
[5; p.225] and (iii) a corresponding disturbance- state 
observer (both full-order and reduced-order) have 
turned out to comprise a surprisingly-effective tool in 
the design of “smart”, high-performance, modem 
MIMO control systems [6]. In particular, such 
control-systems can automatically “accommodate” 
(i.e., cancel-out , minimize or optimally utilize) the 
presence of a broad variety of uncertain, persistent, 
uncontrollable-inputs (disturbance-inputs) w t {t) that 
would otherwise interfere-with the desired closed- 
loop behavior of a multivariable dynamic system, [7]. 

In connection with practical applications of 
disturbance-accommodating control system design 
many users have observed that a disturbance- 
observer is capable of adapting-to, and producing 


surprisingly-close, real-time estimates w(t) of, the 
actual disturbance time-behavior w(t) - - - even when 
the actual time-behavior of w(t) deviates from the 
behavior modeled by the “internal-copy” [5; Fig. 2 
and text below] of the w(t) state-dynamics embodied 
in the disturbance-observer design (i.e., in the spline- 
model and D-matrix). 

In this paper a technical explanation, with 
supporting examples and simulation results, is 
presented to explain the adaptive/robust disturbance- 
estimation capability of a disturbance-observer. 

II. Spline-Models and State-Models for Uncertain 
Disturbance-Inputs 

To minimize the complexity of this presentation, 
we will restrict attention to the simpler case of a 
single (and scalar) disturbance-input w(A However, 
in principle, all our results apply also to the case of 
multivariable disturbance- inputs 

w(t) = w = col.(w I (0,w 2i (0»"-’»Wp(0)- 

In the mathematical-theory of disturbance- 
accommodation, [7], [8], the uncertain time-behavior 
of a (real-valued, scalar) disturbance-input w(t) is 
modeled by a spline-function of the form 

Mt)=C l f l (t)+C 1 f 2 (t)+---+C M f M (t); w= scalar, (1) 

where {c^ is a set of “arbitrary,” real-valued 

weighting-coefficients that are presumed essentially 
“constant” but may abruptly change-values in a time- 
sparse, once-in-a-while manner (hereafter referred-to 
as a “stepwise-constant” manner). In (1) the set 

{/Mr ls a set rea F valued, well-defined, well- 

behaved, (continuously differentiable, etc.) 
independent, known, “basis-functions” that are user- 
chosen to represent the user’s knowledge/ 
presumptions about the fundamental-modes 
(building-blocks) of w(t) time-behavior. At each 
moment the basis-functions are “weighted” and 
linearly-combined according-to (1) to produce the 
actual time-behavior (kinematics) of w(t) at that 
moment. 

In practical applications the time-behavior of w(t) 
is “uncertain” because the values of the “stepwise- 



constant” weighting coefficients in (1) are 

(typically) unknown, not reliably predictable, and not 
directly measurable. Moreover, the user’s a priori 

choice of basis-functions might not exactly 

“match” the actual basis-functions associated with 
real-life, w(t) time-behavior encountered in a 
realistic, practical application. 

In disturbance-accommodation applications, 
the practicality and effectiveness of the spline-model 
(1) is considerably enhanced if the user chooses 
basis-functions f(t) in (1) from the class of 
functions that are of “linear-dynamic” (LD) type, [8; 
Eqs. 6-12]. That is, each fi(t ) obeys a respective 
well-behaved , well-defined linear, homogenous 
differential equation (of finite-order). In that case, on 

those sub-intervals Q of time where the all 

remain constant, the function w(t) defined by (1), 
with “arbitrary” (constant) weighting- 

coefficients {C i , also satisfies a knowable, linear, 

homogenous differential equation (of finite-order) of 
the form 


d p w d p l w dw ■ _ 

+ cc r + ’” + a~ Ha,w = 0, t g Q, 

df p dr 1 2 dt 1 


( 2 ) 


where the in (2) are known/knowable 

parameters that depend only on the basis set {^(0}^ 

and, in some cases, may be time-varying parameters, 
[8]. Consequently, the one, higher-order differential 
equation disturbance-model (2) can be replaced by an 
equivalent “disturbance state-model” [6] consisting 
of a set of “p” first-order , ordinary, coupled 
differential equations of the form 

&t=gi(z (*)=; (3) 


where the “disturbance state-variables” z, are 
suitably-defined independent, functions of w{t) and of 
the first (p-1) time-derivatives of w(t). For instance, 
as one illustration, one can define the z , as the “phase- 
variables” 

z x =w, ; z 2 =w ; z 3 = w, z p = w , (4) 


in which case the disturbance state-model (3) 
becomes (in light of (2), (4)) 

i* = l,2,--,(p-l), 

U=-a 1 z 1 -a 2 z 2 — -a p z p . 


The results (4), (5) can be more-compactly expressed 
in the vector-matrix format 

w(t)=<h o ,z(t)> ; <*,* > -denotes inner product (6a) 
z = D 0 z \ z(t Q ) = arbitrary, teQ (6b) 

where 

h o = row(l,0,0, •••,()) (6c) 
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As shown in [6; p.636] the restricted “disturbance 
state-model” (6) can be easily generalized to 
correctly represent the time-sparse, unpredictable 
jumps in the values of the “constants” C, in (1) by 
adding, to the right side of (6), the p-vector input- 
term o(t) to obtain z = D a z + a(t) where 

ct(0 = co1.(ct,(0,ct 2 (/)•••, <T p (0) , (7) 

and where each element o t (t) of a(/) is a time-sparse 
sequence of uncertain dirac-impulses having 
unknown intensities and unknown (sparse) arrival- 
times. The effect of the impulse-sequences aft} on 
the solutions z(0 of (6) is to impart uncertain, time- 
sparse jumps in the values of the corresponding z,(/j 
and thereby “cause/model” similar, corresponding, 
unknown time-sparse jumps in the values of the C , in 
(1), [8; p.395]. 

Since the primary purpose of this paper is to 
provide a technical explanation for a disturbance 
observer’s apparent, inferred “adaptive-behavior” of 
the Cj -values in (1), between successive time-sparse 
jumps in the C„ (i.e. for t e £1) it suffices to consider 
the corresponding restricted state-model (6) in the 
remaining sections of this paper. Moreover, to 
further simplify things, we will hereafter assume the 

parameters {a, }^ in (2), (6) are all constant (and real- 
valued). 

III. Real-Time Disturbance-State Observers 

It is recalled that the conversion of the LD spline- 
model (1) into a homogenous, linear differential 
equation model (2) and subsequently into a 



homogenous, linear, “disturbance state-model” (6) is 
predicated on the tacit assumption that the set of 

“stepwise-constant” weighting coefficients {cff in 

(1) remain constant over the time-intervals t e Q. 
being considered. That simplifying assumption 
allows us to avoid several irrelevant side-issues that 
would distract from the previously stated primary 
purpose of this paper. 

To develop the results in this paper, we lose no 
generality (and achieve further reduction of 
complexity) by employing the simplified, continuous- 
time, disturbance-state observer introduced-in [8; p. 
434] corresponding to the special (and relatively rare) 
case where the (scalar) disturbance w{t) is directly, 
reliably and continuously measurable by some 
appropriate (noise-free) sensor. For that special case, 
the “full-order” version of a continuous-time, 
“disturbance-state observer”, as presented in [8] and 
based on (6), has the simple form [8; Eqs. (52), (53)]. 

Z = D 0 Z- k d0 [w(t) - (h 0 , £(?))], (8) 

where w(t) in (8) denotes the real-time measurement 
of the actual disturbance (i.e. the “input” to the state- 
observer (8)), and k do denotes the observer gain- 
vector that must be designed to assure 
z(t) -» z(t) sufficiently fast for the application 
requirements. For the latter purpose, the dynamical 
equation governing the observer-error e z = (z - £) is 
easily computed to be [8 ; Eq. (53)] 

4 =\Po +WK » (9) 

and thus the gain-vector k do should be designed to 
make all solutions e z (t) of (9) approach zero 
sufficiently fast. Since we are assuming D a is a 
constant matrix, it suffices to design k do to be a 
constant vector that makes all eigenvalues of the 
composite-matrix \D 0 + k do h 0 ] lie sufficiently-deep 
in the left-half of the complex-plane. 

IV. Relation Between the {C'}^ in (1) and the 
{Zj}[ in (6) 

In the routine practical applications of 
disturbance-accommodating control-theory, it is 
neither necessary nor beneficial to “know” the values 
of the Q in (1), thanks to the ability to express the 
best/optimal control-input (control-law) in terms of 
reliable, real-time observer- estimates of %(t) (and of 
plant states x^t)). This fortunate “ability” is so 
important it is called the Principle of (Real-Time) 
Optimal Disturbance- Accommodation [7; pg. 219- 
20 ] 


However, since the spline-model (1) and the state- 
model (6) both represent the same (scalar) 
disturbance time-behavior w{t), it is clear that the 
values of the Cj in (1) must be related-to the values 
of the state- variables z; in (6). In fact, if the values of 
the {Zj{t)f x , p < M, are known at any moment t e Q, 
then (in principle) the corresponding values of the 
in (1), at that same moment, should be 

uniquely determined (recall the basis-functions fi(f) in 
(1) are assumed independent). In the LD case, as 
considered here, the Q are related to the Zj by a set of 
coupled, linear algebraic equations; see [9; Eq. (5d)] 
and the Examples below. 

This observation suggests a novel procedure for 
experimentally evaluating the real-time values of the 
Q in (1) from reliable, real-time state-estimates 
Zj(t) obtained from a disturbance-state observer 

if and when there is some “need” to know/monitor 
the real-time C r values in (1). 

V. What Happens if the Actual w(/)-Behavior 
Doesn’t Conform to the Models (1), (6)? 

Many users of disturbance-accommodating 
control theory have “observed” that disturbance-state 
observers (both full-order and reduced-order; 
continuous-time type) seem to correctly estimate the 
actual disturbance-input w(t) even when the actual 
w(t) time-behavior (mildly) deviates from the 
modeled time-behavior embodied in the spline-model 
(1) and in the state-model (6). For instance, if (1), (6) 
are based-on the behavior-model: w(t) = unknown 
constant = C u the corresponding observer (8) (with a 
sufficiently short “settling-time) produces real-time 
estimates w(t) = h Q z(t) that are surprisingly accurate 
even when the actual w(t) is not constant but is 
(slowly) varying with time. The traditional intuitive 
explanation for this seemingly “adaptive/robust- 
estimation behavior” is that a continuous-time 
disturbance-observer is continually re-evaluating, the 
state-estimate z(t), and reconciling z(/),with (1) 

based-on “current” actual w(0-behavior and on the 
models (1), (6). Thus, in the case cited, a slowly- 
varying actual w(t) is “seen/perceived” by a 
continuous-time disturbance-state observer as a series 
of (infinitesimally) short, “constant” stair-step 
variations \ in the expected w(f) = constant behavior 
embodied in (1), (6) and in the corresponding 
observer D-matrix structure. 

In the remaining sections of this paper we will 
experimentally “confirm” this intuition by means of 
several worked examples using, realistic “composite- 
state” (plant + disturbance) observers [8; pp. 431- 
437], and corresponding numerical, simulation- 
results; In particular, we will use the algebraic 



equations relating the Q in (1) to the z f in (4)-(6) to 
experimentally evaluate the disturbance observer’s 
“perception” (re-construction) of the actual, real-time 
values of the “weighting-constants” C t in (1) when 
the actual w(t) time-behavior (re Q) deviates (both 
“mildly” and “grossly”!) from the time-behavior 
modeled by (1) and embodied in the structure of the 
disturbance observer’s D-matrix in (8). These results 
demonstrate that when the actual w(f)-behavior 
deviates from the behavior modeled by (1), the state- 
estimate z(t), generated by a continuous-time 
“disturbance-state observer” infers that the 
“constants” Q in (1) are rapidly changing-value in a 
virtually-continuous, stair-step manner, so that the 
resulting w(t) in (1) then closely “matches” the actual 
w(0-behavior. 


performance of the disturbance observer when an 
appropriate dynamic basis for the observer is chosen 
when a “poor” choice is made for the dynamic basis 
of that observer. Let us first consider an example 
where the basis chosen for the disturbance observer is 
a “good” one (“good” in the sense that it contains a 
basis representative of the disturbance). 


Example #1: 

• Disturbance a constant with a step-change at 
t = 3 s. 

• Spline-model for the disturbance observer is 
w(t) = cj + c 2 t (a “good” model) 

The disturbance state model is: 




0 1 
0 0 


H=[l,0] 


(11) 


VI. Illustrative Examples 
To demonstrate the apparent “adaptive-behavior” 
exhibited by disturbance-observers a simple pure 
inertia plant will be modeled with a disturbance 
applied. The system layout with a composite 
observer is given in Fig. 1. 
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Figure 1: Overview of the Simulation 
The state-space representation of the plant model is: 
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cancels w(t) effects. 


State feedback control u s - K p x is achieved by 
multiplying the state estimate X by the gain Kp. 
Disturbance cancellation is achieved by feeding back 
the disturbance estimate w to achieve B p u d ~-Fw. 
The the closed loop roots of the system at are placed 
at [-2 -2.5]. This results in Kp = [5 4.5] This value 
of Kp will be used throughout the examples presented 
below. 

To demonstrate the “adaptive” behavior of the 
disturbance observer we will first examine the 


Placing the composite observer poles 
[8; p.4-31] at -6,6±2.4i and -6.9±1.2i results in the 
following composite observer feedback gain: 

Ko = col(26.9, 279, 1 320, 2400) (12) 

The results of this example are shown in Fig. 2. 


Observer Basis: c 1 + c 2 1 



Fig. 2: Constant Disturbance, Cj + c 2 t Basis 
The disturbance w(t) has an initial value of -0.5 and 
at 3 sec jumps to a value of 1.0. As one might expect 
(since the observer dynamic basis is well chosen for 
this disturbance) the composite observer does a good 
job of estimating the disturbance and the observer 
identifies the c h c 2 as “stepwise-constants”. When an 
unpredictable jump occurs at t = 3 sec the observer 
experiences a transient but quickly recovers and locks 
back onto the disturbance w(t). In Fig. 2 (a) we see 
the disturbance state z(t). Since H = [1 0] the 
disturbance estimate w(t) = £, (/) (6a). Examination 
of zj(t) in Fig. 2 (a) shows the good tracking of the 




disturbance w(t) spoken of above. Additional insight 
can be gained by looking at the observer performance 
in terms of weighting coefficients in (1) for the spline 
function w(t) = Cj + c 2 t. It is simple to solve for c } 
and c 2 in this example (this is not always so easy): 


Cy = z { - z 2 t 


(13) 


C 2 = Z 2 

Fig. 2(b) shows that the c 2 coefficient has a 
significant role during the transient and then returns 
to zero as the observer locks back onto the constant 
disturbance and only Cj is needed. Fig. 2 (c) confirms 
that the spline model of the disturbance 


w(t) = Cj (t) + c 2 (t)t is identical to the 


disturbance-state representation w(/) = z x {t) . 


Example #2 

• Disturbance is a sinusoid with a frequency 
of 0.75 rad/sec and an amplitude of 1 .0. 

• Basis for the disturbance observer is 
W(t) ~ Cj + c 2 t 

The disturbance state model is the same as in 
example #1 (11). The observer poles are also at the 
same location so the observer gain Ko is the same 
(12). This observer has the same basis so the relation 
between z(t) and c(t) is also the same (13). 

Observer Basis: c 1 + c 2 1 



Fig. 3: Sinusoid Disturbance, Cj + c 2 t Basis 


• Spline-Model for the disturbance observer is 
w(t) = cj sin(cot) + c 2 cos(cot 


The disturbance state model defined in (6) for this 
disturbance basis is: 



1 

0 


H=[o> 0] 


(14) 


Composite observer poles at -6. 6±2.4i and - 
6.9±1.2i result in an observer feedback gain of: 

Ko = col(26.9, 279, 1740, 2990). (15) 

The relationships between the disturbance-state 
variables z(t) and the spline coefficients c(t) for this 
observer are: 

c 2 = coZj cos(co0 - z 2 sin(coO 
^ = (cozj - c 2 cos(co/)) / sin(coO- 


In this example we are tracking a constant 
disturbance (with a jump) using an observer designed 
for the ; spline-model Cj sin(cot) + c 2 cos(cot). One 
would expect poor disturbance tracking but the 
results shown in Fig. 4 demonstrate surprisingly good 
tracking performance. 



Fig. 4: Constant Disturbance, Sinusoidal 
Basis 


The surprising result shown in Fig. 3 is that a 
disturbance observer designed with a dynamic model 
(“internal-copy”) for a constant plus a ramp is able to 
track a sinusoid quite well. Both c } and c 2 
continually-adapt their “constant-values” to maintain 
accurate tracking of the disturbance. 

Example #3: 

• Disturbance a constant with a step at t = 3 s. 


Note that in Fig. 4 the values of and c 2 are 
constantly shifting value/”adapting” as the sine and 
cosine rise and fall with time “t”. The behavior can 
be explained by the fact that the disturbance observer 
constantly attempts to minimize the error between the 
disturbance estimate w(t) and the actual disturbance 

w(t) (8). Even if the basis of the observer dynamics 
in no way reflects the dynamics of the actual 






disturbance the inherent property of a state-observer 
is that it will, according to the error the dynamics in 
(9), seek to minimize the error between the 
disturbance and its estimate. This property of an 
observer gives a control system gives what can be 
portrayed as “disturbance-adaptive” behavior. 

There are limitations to how well an observer can 
track a disturbance with a “poor” selection of the 
observer basis dynamics. The performance of a 
sinusoid based disturbance-observer deteriorates as 
the frequency © used in the disturbance-state model 
(14) becomes large relative to the frequency of jumps 
in the “constant” type disturbance. Similarly, the 
performance using a constant plus ramp basis for an 
observer tracking a sinusoid deteriorates as the 
frequency of the sinusoidal disturbance increases. 
This can be offset somewhat by placing the observer 
poles deeper into the left-half-plane. 

VI. Conclusions 

The tracking performance of disturbance-observers 
when the internal-copy of w(fy in the observer is not a 
good representation of a given disturbance dynamics 
has been studied. The ability of a disturbance- 
observer to track a disturbance with dynamics 
completely different than what was assumed by the 
observer designer is the somewhat surprising result 
that this paper documents. The beauty of using a 
disturbance-observer in the control scheme is that 
even with totally unanticipated disturbance dynamics 
the “adaptive” capabilities inherent in the observer 
can still yield good results. This observer “adaptive” 
feature is not universal (the control designer should 
always try to include all important disturbance 
dynamics in the disturbance-observer basis) but it 
does account for some of the robustness that 
characterizes control systems that use a disturbance- 
observer to accommodate disturbance inputs. 

Acknowledgement: Thanks to Dr. C.D, Johnson for 
his assistance in the preparation of this paper. 
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Introduction 


• A disturbance observer is a mathematical tool used to 
estimate in real-time the state of a disturbance acting on 
the system. 

• With knowledge of the disturbance state a controller can 
compensate (accommodate) for the disturbance. 

• We know that when the math model for the anticipated 
disturbance is a “good” one that a disturbance observer 
can quickly acquire and track the disturbance. 

• It has been observed, however, that these observers 
display unexpected robustness when the functions 
chosen as the basis set for the observer are a poor 
reflection of the disturbance waveform. This robustness 
is what this presentation will illuminate and try to explain. 




Background 











Background 


• Many disturbances have waveform structure with once-in-a-while 
jumps 

- As a the discontinuous jumps in a disturbance occur faster & faster it 
approaches “noise” and statistical methods should be used. 

- Many real world disturbances have waveform structure with only 
occasional jumps and can readily be estimated with a disturbance state 
observer 

• A composite observer estimates both the plant state and the 
disturbance state 

- To estimate the disturbance state we must first decide on an 
appropriate disturbance model. Common disturbance state models 
include: 

• Constant 

• Constant + ramp 

• Sinusoid 

• Pulse 

• Some combination of the above 
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Waveform Model of a Disturbance 


• We must first mathematically model the anticipated 
disturbance “w” 

w (0 = + ••• + c M f M (t ) (1) 

• The f t (t) are known independent functions of time 

• The c i are unknown weighting coefficients 

• This is referred to as the function space representation of w(t) 

• This set of functions f/t) form the basis-set and is 
chosen to reflect the expected waveform patterns of the 
disturbance 

• The waveform model (1) is the mathematical foundation 
for this approach to disturbance modeling 

- The waveform representation, however, can not be 
readily implemented into a control algorithm 




Disturbance State Models 


• In order to use a disturbance observer one must first go to an alternative 
representation 

- This alternative to a waveform representation is a disturbance state model 

• Determining a disturbance state model can be approached as a reverse 
differential equations problem. The waveform representation in (1) can 
also be represented as the solution to the following general DE 


d p w 

dt p 


+ cc . 


d p ~'w 

dt p 


+ • • • + a 2 


dw 

dt 


+ a x w - 0 



- The (Xj in (2) are known parameters that depend only on the basis set f/t) 
• Equation (2) can be written in the vector-matrix format: 


w{t) = Hz 



z = Dz + cr(t) 


The <r(t) term represents a series of dirac-impulses having unknown intensity and unknown arrival time. 
This unknowable term is not part of the design but is included for mathematical completeness. 
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Example Disturbance Model Representation 


• For example, let us construct a disturbance state model for a 
disturbance that is mostly constant with occasional jumps 



A disturbance w(t) that is constant 
with “occasional” jumps 


Waveform Representation: w(t) = c x 


Disturbance State Representation: 


dw 

~dt 



H = [ 1], D = [0] 
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Example Disturbance Model Representation 


• For a periodic disturbance with changing amplitude but a fixed 
frequency calls for a sinusoid basis function 

A disturbance w(t) that is periodic 
with changes in amplitude 


Waveform Representation: w(t ) = c, sin cot + c 2 cos cot 

1 2 

Disturbance State Representation: — — + a, w= 0 

dt 2 1 



H = [6> 0], Z) 


0 1 

-co 1 0 
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Simulation Example 


Will use the control of a pure inertia system as our test-bed 

Will examine the performance of the disturbance observer for 
“good” and “bad” choices of basis functions 


Constant 


u = u p + u d 
u p regulates x(t) 
u d cancels w(t) 


Disturb 

estim 


S 








~0 

r 


'O' 


'o' 

A = 

0 

0 

B = 

1 

F = 

1 







~ 1 


c = [1 0 ] 



Plant Model 
dx = Ax + Bu + Fw 


x vec 



xHat 

2 


x estim 
2 





output 


wHat 


Composite Observer 

dX est = AX est + K °<y- Cx est> + Bu 


w estim 
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Example 1 


Disturbance is a constant with an unexpected jump 
Disturbance observer: “constant + ramp” basis (a “good” choice) 


Disturbance Model : w(t) = c, + c 2 t D = 


0 1 
0 0 


H = [1 0 ] 


Observer Basis: c 1 + c 2 1 




Time (sec) 


With a constant in the 
disturbance observer basis set 
one would expect the good 
tracking seen in the Figure 
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Example 2 


Disturbance is a sinusoid (oo = 0.75 rad/s) 

Disturbance observer: “constant + ramp” basis (a “poor” choice) 


Disturbance Model : w(t ) = c, + c 2 J 


D = 


Observer Basis: c 1 + c 2 1 


0 1 
0 0 


H =[10] 




The observer tracks the 
sinusoid disturbance using a 
constant + ramp internal model 
for the disturbance. Note how 
c 1 & c 2 constantly adjust to 
minimize the error between the 
disturbance and its estimate 
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Example 3 

• Disturbance is a constant with an unexpected jump 

• Disturbance observer: Sinusoid basis (a “poor” choice) 


Disturbance Model : w(t) = c x sin cot + c 2 cos cot 
Observer Basis: c^intcot) + c 2 cos(cot) 



1 

0 


H = [co 0] 





An observer designed for a 
sinusoid disturbance successfully 
tracks a constant disturbance! 

The coeffs c 1 & c 2 are rapidly 
changing but the net disturbance 
estimate (14;) smoothly & accurately 
tracks the disturbance 


(C) 



Time (sec) 
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Conclusions 


• Disturbance observers demonstrate an ability to “adapt” and 
successfully track a disturbance with dynamics completely different 
from what the observer was designed for 

- There are limits to this adaptability. It is always best to pick a basis-set 
that is representative of the expected disturbance dynamics 

• The ability of an observer designed for a sinusoid to track a constant with 
jumps deteriorates as the frequency of the sinusoid basis becomes large 
relative to the frequency of the jumps 

• Similarly, the performance of a disturbance observer with a constant basis 
function deteriorates as the frequency of the sinusoid disturbance increases 

• An observer is a closed loop subsystem that is structured to drive an 
estimate to zero. As a result of this disturbance observers will make 
every effort to adjust internal coefficients so that the net error is 
minimized. This effect manifests itself as the “adaptive” capability of 
a disturbance observer to track disturbances very different from 
what it was originally designed for. 
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Backup 
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Derivation of the function-space Coefficients 


• The disturbance-state z is readily available from the 
observer. Much insight can be gained by examining the 
dynamics in the function space. 

- The disturbance state z can be converted to the c, values in (1) 
needed to represent w(t) in the wavefrom model 

• Disturbance observer with a “constant + ramp” basis 

Zj = w = q + c 2 t and z 2 = z l = ~ (c, + c 2 t ) = c 2 


Therefore 


C\ “l ^2^ 
c 2 = z 2 
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Derivation of the function-space Coefficients 


Disturbance observer with a “sinusoid” basis 

- Recall the waveform model is w(t) = c, sin at + c 2 cos cot 

- The state model is: 

w(t) = Hz, z = Dz => H = \a> 0 ], D - 


0 


-co 


1 

0 


Start with: w(t) = coz l = c x sin cot + c 2 cos cot 

c c 

z x = — sin cot + —cos cot 

CO CO 


Take the derivative of z, — > — z 2 — c, cos cot c 2 sin 

Have 2 equations & two unknowns. Solving we find that 

c 2 = cozj cos(co/) - z 2 sin(co^) 
e l = (oo z x - c 2 cos(co t)) / sin(co/). 
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Example B1 


Disturbance is a sinusoid 

Disturbance observer: Sinusoid basis (a “good” choice) 


Disturbance Model : w(t) = c, sin cot + c 2 cos cot 


D = 




0 1 

2 


-CO 


0 


H = [cd 0 ] 


The disturbance is a pure sine with no 
phase angle. As one might expect the 
c 1 coeff goes to 1 .0 and the c 2 coeff 
goes to zero after the initial startup 
transient. There is a startup transient 
because the initial conditions of the 
observer integrators are set to zero. 
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